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ABSTRACT 

We present mid-infrared spectral mapping observations of the core-collapse supernova remnant 
IE 0102.2-7219 in the SmaU Magellanic Cloud using the InfraRed Spectrograph (IRS) on the Spitzer 
Space Telescope. The remnant shows emission from fine structure transitions of neon and oxygen as 
well as continuum emission from dust. Comparison of the mid-lR dust emission with observations at 
x-ray, radio and optical wavelengths shows that the dust is associated with the supernova ejecta and 
is thus newly formed in the remnant. The spectrum of the newly formed dust is well reproduced by 
a model that includes 3 x 10"'^ Mq of amorphous carbon dust at 70 K and 2 x 10~^ Mq of Mg2Si04 
(forsterite) at 145 K. Our observations place a lower limit on the amount of dust in the remnant since 
we are not sensitive to the cold dust in the unshocked ejecta. We compare our results to observations of 
other core-collapse supernovae and remnants, particularly Cas A where very similar spectral mapping 
observations have been carried out. We observe a factor of ^ 10 less dust in E0102 than seen in Cas 
A, although the amounts of amorphous carbon and forsterite are comparable. 
Subject headings: dust, extinction — infrared: ISM — supernova remnants 



1. INTRODUCTION 

Dust is a crucial component of the interstellar medium 
(ISM). It provides a site for interstellar chemistry, in par- 
ticular the formation of molecular hydrogen; it regulates 
thermal balance in various phases of the ISM; and it pro- 
vides shielding for dense clouds, critically influencing the 
process of star-formation. The characteristics of inter- 
stellar dust depend on the balance between its formation 
and destruction and on the processing that dust grains 
endure in the ISM. The mechanisms for dust production, 
their tiniescales and efficiencies and the effect of metal- 
licity on these quantities determine the cosmic history of 
dust production and are mostly not well characterized. 

In particular, the relative amount of dust that core- 
collapse supernovae (CCSN) contribute to the ISM is 
poorly constrained and a matter of much debate. Be- 
cause of the short lifetimes of their massive progenitor 
stars (a few tens of Myrs), CCSN are often cited as an im- 
portant source of dus t early in the histo ry of the Universe 
ilMorgan et al.ll2003l : lOwek et al.l[2007HTodini fc Ferraral 
|2P0l). Observational evidence for substantial dust pro- 
duction in supernovae, however, is scant and controver- 
sial. 

Our understanding of dust formation in the early Uni- 
verse and the role of CCSN has recently been challenged 
by observations of large dust masses around high rcdshift 
quasars (Bertoldi et al. 2003; Bcel en et all [200 6). Be- 
cause the redshifts of these quasars are so high and the 
dust masses are so large, the dust had to be produced 
efficiently and on a short timescale. For the z — 6.42 
quasar SDSS J1148-I-5251, the age of the Universe is 
only ~ 900 Myr and more than 10^ Mq of dust are 
observed in its vicinity (although there are substantial 
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uncertainties associated with this dust mass, see the dis- 
cussion in Omont ct al. (2001)). There is some indication 
that dust in some of these high redshift galaxies may fol- 
low an extinction curve unlike th at in either the Milk y 
Way or Small Magella nic Cloud (jMaiohno et al.ll2004D . 
iHirashita et al.l ()2005l ) claim that the extinction curve 
of dust in the z = 6.19 quasar SDSS J1048+4637, for 
example, can be reproduced by the species and size dis- 
tribution of dust formed in models of the unmixed ejecta 
of ^ 10 — 30 Mq supernovae. 

The models of dust production in supernovae typi- 
cally predict on the order of 0.1 — 1 Mq of dust formed 
in an average CCSN (see fo r instance, [T odini fc F erraral 
120011: iNozawa et al.l [20031) . iDwek et alT(r2007 ) estimate 
that to produce the mass of dust observed around SDSS 
J 1148-1-5251 requires an average yield of ~ 1 Mq of 
dust per CCSN. As of yet, however, observations of 
newly-formed dust in supernovae and their remnants 
have shown orders of magnitude less dust than would be 
expected based on these models. If we cannot substan- 
tiate the predicted dust production efficiency of CCSN, 
our explanation of the high-z quasar observations may 
need to rely more heavily on the contributions of red 
supergiants and mass i ve AG B stars to dust production 
(|van Loon et al.|[2005l [200l . 

Part of the discrepancy between expectations and ob- 
servations of dust production in CCSN and their rem- 
nants stems from the difficulty of observing newly-formed 
dust. The condensation of dust is thought to happe n be- 
tween ^ 400 — 800 days after the explo sion (Kozasa et al.l 
[T989l:lTodini fc Ferraral200l[ : INozawa" et al...200a) . While 
the newly-formed dust is still warm, its near- and mid- 
IR emission may be detectable. Observations of the 
thermal emission newly-formed from dust in SN 1987A 
show ^ 10^"* M0 of dust produced by 775 days after 
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the explosion (iRoche et all 119931: IWooden et al.l 119931 : 
iMoselev et all 119891 ) . ISugerman et all (|2006f ) claimed 
to detect mid-IR emission from ~ 0.02 Mq of newly- 
formed dust in_SN_2003gd, althou gh this claim has been 
disputed by 'Meik le et alJ ([2007) after their reanalysis 
of the Spitzer observations that showed only 4 x 10^^ 
Mq of dust. At these early times, near- and mid-IR 
emission can also arise from IR light echoes, as the su- 
ernova's radiation heats pre-exis t ing circumst ellar d ust 
Bode fc Evai^[l980l : iDweafTosl [Po^o et al.ll2004D . If 
there is a contribution by an IR light echo, the mass of 
newly-formed dust inferred from the early-time infrared 
emission may be greatly overestimated. On the other 
hand, dumpiness of the ejecta could lead to an underes- 
timate of the dust m ass based on the infrared emission 
(|Ercolano et al.ll2007| ). 

As the supernova ejecta expand and cool, the newly- 
formed dust becomes detectable only at far-IR and/or 
submillimeter wavelengths at which time confusion with 
cold dust in nearby sta r- forming regions becomes a ma- 
jor issue. iDunne et al.l f2003.) claimed a detection of 2 
M0 of newly- formed dust in Gas A with 850 /xm SCUBA 
observations. iKrause et al.l ()2004f ). however, used far-IR 
and molecular line observations to argue that there is 
at least an order of magnitude less newly-formed dust 
in Cas A and that most of the submillimeter emission 
arises in an intervening molecular cloud. Similarly, sub- 
millimeter observatio ns of Kepler's remnant were used by 
iMorgan et all (|2003f ) to argue for 0.1 - 3 of newly- 
formed dust although these results have been recently 
disputed by Blair ct al. ( 2007) who find evidence for only 
5 X 10"'' Mq using far-IR observations from Spitzer. 

Eventually, the interaction of the ejecta with the 
surrounding circumstellar and/or interstellar medium 
(CSM/ISM) produces a reverse shock which propagates 
back into the ejecta. The reverse shock reheats the 
newly-formed dust so that once again it is detectable in 
the mid-IR and may destroy some of the dust via sputter- 
ing, as well. Observations at 24 ^m of IE 0102.2-7219, 
the sa me SNR we discuss here, by IStanimirovic et aTl 
(|2005D showed evidence for - 8 x lO"'' Mq of dust. Re- 
cently, iRho et aTl ()2008f ) used IRS on Spitzer to make a 
spectral map of Cas A. They claim evidence for ^ 0.02 
Mq of newly- formed dust in the remnant. Since only 
the reverse-shocked dust should be visible in the mid- 
IR, only a fraction of the newly-formed dust can be de- 
tected in this way. Additionally, circumstellar or inter- 
stellar dust that is interacting with the forward shock 
should also be detected in the mid-IR. I Williams et al.l 
(|2006l ) observed four supernova remnants in the Mag- 
ellanic Clouds with the IRAC and MIPS instruments 
on Spitzer and concluded that the mid-infrared emission 
was from shocked circumstellar dust rather than newly- 
formed dust. Despite these limitations mid-IR observa- 
tions are very useful for studying the dust production 
in supernovae because they typically have higher an- 
gular resolution, which allows the remnant to be more 
easily separated from foreground and background emis- 
sion. Mid-IR spectroscopy also has the potential to pro- 
vide some information on the composition of the newly- 
formed dust in a remnant from the distinctive mid-IR 
spectral features of some dust species. In this paper 
we present mid-IR spectral mapping observations of the 
Small Magellanic Cloud remnant IE 0102.2-7219, which 



we use to learn about the mass and composition of dust 
produced in CCSN. 

SNR IE 0102.2-7219 (E0102) is a young, oxygen-rich 
remnant located on the outskirts of the emission-line neb- 
ula N 76. It is thought t o be the product of a Type Ib/Ic 
or IIL/b supernova ( Blair et all l200fll : iChevaheii l2005h 
due to the lack of hydrogen and helium in the remnant's 
spectrum and the high abundance of oxygen and neon 
in its ejecta. E 0102 has been observed extensively at x- 
ray, UV, optical, infrared and radio wavelengths. Chan- 
dra x-ray imaging of the remnant by iGaetz et al.l ([2000,) 
showed a textbook remnant structure with a blast wave 
radius of ^ 22" an d a reverse shock radius of ~ 15" . 
lAmv fc: Balll (jl993f ) observed the remnant at 6 cm us- 
ing the Australia Telescope Compact Array (ATCA) and 
found a shell of radio emission with a radius of ~ 20". 
This ring lies between the measured radii of the faint rim 
of the blast wave and the bright ring of the reverse shock, 
and is attributed to material in the forward shock. The 
full-width velocity dispersion of the optical [O III] emis- 
sio n from the remnant was measured to be ~ 6000 km/s 
bv iTuohv fc Dopital (|1983D . im plying an age of ~ 1000 
years given the remnant's size. iHughes et al.l ()2000f ) use 
three epochs of x-ray observations from the Einstein, 
ROSAT and Chandra observatories to measure an ex- 
pansion rate for the whole remnant in agreement with 
the results of iTuohv fc Dopital (|1983l ). Recent work by 
iFinkelstein et al.l l 20061 ) used two epochs of optical [O 
III] observations to determine proper motions for ejecta 
filaments and find an age of 2050 ± 600 years for the rem- 
nant. The results of the following work do not depend 
sensitively on the age, so we adopt a value of 1000 years 
f or simplicity. 

IStanimirovic et aTl (|2005( ) presented the first infrared 
detection of SNR E 0102 from the Spitzer Survey 
of the Small Magellanic Cloud (S^MC) observations. 
S^MC covered the star-forming regions of the SMC Bar 
and Wing with ima ging in all MIPS and IRAC bands 
(jBolatto et al.ll2007f). The remnant was on ly clearly de- 
tected at 24 fim. IStanimirovic et all (2005) used the 24 
fim detection along with the upper limits in the other 
bands to argue that the dust in the remnant has a tem- 
perature of ^ 120 K and they determined an upper limit 
of 8 X lO^'' Mq of warm dust. The main uncertainties 
in their analysis were the fraction of the MIPS 24 /im 
emission coming from the [O IV] line at 24.8 /im and the 
temperature of the dust. 

In the following, we present Spitzer IRS spectral map- 
ping observations of SNR IE 0102.2-7219 obtained as 
part of the Spitzer Spectroscopic Survey of the Small 
Magellanic Cloud (S'*MC). These observations allow us 
to map the emission lines and dust continuum from the 
remnant and place a lower limit on the amount of newly 
formed dust in E 0102. In Section l2Tl we describe the ob- 
servations and their reduction. In Section [3. II we present 
the spectrum of the remnant and discuss the spatial dis- 
tribution of its various components. In Section 13.31 we 
use a model for the reverse-shocked dust to measure a 
mass of 3 X 10"'^ Mq and we discuss the composition of 
the dust as revealed by its spectrum. Section |4] discusses 
the implications of our measurement for the understand- 
ing of CCSN dust production and compares our results 
with previous work. 
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2. OBSERVATIONS AND DATA REDUCTION 
2.1. Spitzer Observations 

SNR E 0102 was mapped with all LL and SL orders 
of IRS on the Spitzer Space Telescope as part of the 
S^'MC project (Spitzer Spectroscopic Survey of the Small 
Magellanic Cloud, GO 30491). The wavelength coverage 
of the SL and LL orders extends from 5.2 to 38.0 fj,m 
with a spectral resolution between 60 and 120. Ramp 
times of 14 and 30 seconds for SL and LL were chosen as 
a compromise between sensitivity, spatial coverage and 
observation length. 

SNR E 0102 was covered in the spectral maps of the 
N 76 region, observed on 9 and 12 December 2006 for 
SL and LL, respectively. The LL map consists of 6 x 75 
pointings covering an area of 474 x 381 arcseconds^ in 
both LLl and LL2. The SL observations consisted of 
5 X 120 pointings covering 260 x 222 arcseconds^ for the 
SLl and SL2 orders. The roll angle was left uncon- 
strained for our observations. The SL and LL maps have 
different orientations because of the slit positions in the 
focal plane of the telescope. Figure [T] shows a three-color 
image of the region from the S'^MC dataset overlayed 
with the coverage of the SL2, SLl, LL2 and LLl maps. 
The spectral maps are completely sampled in both LL 
and SL by stepping perpendicular to the slit by one-half 
slit width (5.08 and 1.85" for LL and SL, respectively). 
In the SL orders we have stepped by a full slit width 
parallel to the slit (52") in order to increase the spatial 
coverage of the maps at the expense of some redundancy. 
In the LL orders we have used half-slit width steps par- 
allel to the slit (79"). More redundancy in the LL maps 
is necessary due to the increasingly large numbers of bad 
pixels at long wavelengths. 

2.2. Data Reduction 

The spectral maps were calib rated and assemb led in 
IDL using the Cubism package ([Smith et al.l 120071 ). Cu- 
bism takes as input pipeline-processed two-dimensional 
spectral images from IRS and produces three dimensional 
spectral data cubes using a reprojection algorithm which 
employs two steps of polygon clipping (in which the exact 
geometrical overlap of input pixels, or pseudo-rectangles 
in the case of the IRS spectral images, is used to calculate 
its weight in an output pixel grid). Cubism also applies 
an aperture loss correction function (ALCF) and a slit 
loss correction function (SLCF) to adjust the calibration 
which is based on point-source observations rather than 
extended sources. The output pixel scales of the LL and 
SL cubes are 5.08" and 1.85", respectively. Images of the 
IRS data in this work are shown on the IRS pixel scale 
rather than aligned to R.A. and Dec to avoid unnecessary 
interpolation. 

The Galactic and zodiacal light foregrounds were re- 
moved from the IRS spectra using a dedicated "off" po- 
sition. The same position was used for both the LL and 
SL modules. The background observations were done 
immediately following the mapping observations. The 
"off" position was chosen as a place in the S^MC 24 /xm 
MIPS mosaic with no extended emission or point sources, 
located at R.A. l'>09"H0!00, Dec. -73°31'30'.'00. 

The subtraction of the "off" position, in addition to re- 
moving the Galactic and zodiacal foregrounds, mitigates 
the effects of rogue pixels (i.e. time variable bad pix- 



els) in the IRS detectors, although it does not eliminate 
them. Further bad pixel removal was done in Cubism by 
examining each wavelength for striping (caused by global 
bad pixels) and other artifacts. The effect of bad pixels 
increases at the longer wavelengths, degrading the data 
quality beyond about 35 fim. 

2.3. Resolution Matching and Alignment of Cubes 

In order to directly compare the spatial distribution of 
supernova emission at each wavelength as well as prop- 
erly extract a spectrum when the extraction region is on 
the order of a resolution element, it was necessary to con- 
volve all the individual slices of the spectral cube to the 
same resolution. This was achieved by creating convolu- 
tion kernels based on theoretical point-spread functions 
(PSFs) for the IRS spectrograph generated by the pro- 
gram sTinyTim ^ . 

The convolution kernels k were defined in the following 
way: 



FFT^^i?//), 



(1) 



where R is the 2-dimensional Fast Fourier Transform 
(EFT) of the PSF of the longest wavelength shce of the 
LLl cube and I is the EFT of the PSF of the shce in ques- 
tion. The convolution kernels were masked in Fourier 
space such that there was no power at spatial frequencies 
higher than the Nyquist frequency of the lowest resolu- 
tion kernel in order to avoid the amplification of high- 
frequency noise. Each slice of the spectral data cube was 
then convolved with the appropriate kernel. Although 
the PSF of IRS is not well explored (i.e. the predictions 
from sTinyTim are mostly untested) and the PSF that 
Cubism recovers can be somewhat elongated in a way 
that varies with wavelength, we found that the results 
of this technique were adequate for our purposes. The 
convolution process was tested on the cube of the N 22 
region of the SMC which contains a very bright point 
source for which the first two diffraction rings are visible 
for essentially every slice in the cube from SL2 through 
LLl. The correspondence between wavelength slices af- 
ter our convolution is excellent as shown in Figure [2l 
For each cube, the associated uncertainty cube was also 
convolved and the errors appropriately propagated. 

The LLl, LL2, SLl and SL2 cubes were ahgned us- 
ing a technique based on the polygon-clipping algorithm 
implemented in Cubism. In polygon-clipping the weight 
assigned to each pixel in the original image is calculated 
based on the geometric overlap of that pixel with the grid 
of output pixels. Aligning the cubes using polygon clip- 
ping minimizes interpolation error and allows for direct 
propagation of the associated uncertainty cubes. The 
LL2, SLl and SL2 cubes were aligned to match the as- 
trometry and pixel scale of the LLl map after resolution 
matching. 

2.4. Local Background Removal 

More difficult to remove than the Galactic foreground 
is the variable and complex local foreground and back- 
ground emission (for simplicity we will refer to this as a 
background, though there are contributions from the en- 
tire line of sight through the SMC). E 0102 is located near 
the large emission-line complex N76. Emission from the 

^ |http : //ssc ■ spitzer ■ caltech . edu/archanaly/contributed/stlnytlm/inde: 
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Fig. 1. — Three-color image of the N 76 region overlayed with the coverage of the IRS spectral mapping observations. The S^MC MIPS 
and IRAC maps at 24, 8.0 and 3.6 fim are shown in red, green and blue. The white solid and dashed lines show the coverage of the LLl 
and LL2 spectral maps. The cyan solid and dashed lines show the coverage of the SLl and SL2 maps. The IE 0102.2— 7219 supernova 
remnant is located at R.A. 0l'^04™02!l, Dec. -72°01'52'.'5 (J2000). 



outskirts of this region must be removed from the super- 
nova remnant spectrum to accurately measure the dust 
content and emission hne flux. 

An accurate background subtraction is not possible 
over a large area due to the proximity of N 76 and the 
complexity of the emission in the region. For this rea- 
son we focus on the area near the supernova remnant 
and attempt to subtract a model for the local emission. 
We first attempted to fit planar or higher order sur- 
faces to the emission in the region of the remnant. The 
steep profile on the side of N 76 nearest the remnant was 
poorly represented in these fits, resulting in large residu- 
als. Polynomial fits to ea ch row and column, as done by 
iStanimirovic et ahl ()2005D . were a better representation 
of the local emission but also suffered from issues with 
reproducing the steepness of the N 76 shoulder, interfer- 
ing with the goal of determining the background in the 
immediate vicinity of the remnant. Rather than trying 
to reproduce the local emission we argue that a better 
representation of the background comes from interpolat- 
ing across the supernova remnant, which eliminates the 
need to reproduce the complex emission structures with 
polynomials or surfaces. 

For each slice of the convolved and aligned cubes, we 
mask out the immediate region of the remnant, using 



a mask 60" in diameter to exclude any emission from 
the remnant in our background determination. We then 
linearly interpolate in the x and y directions across the 
masked region, using two pixels on either side of the rem- 
nant in the given row or column. Finally, we average 
together the results obtained by interpolating in the x 
and y direction. This technique has advantages and dis- 
advantages. Since we only use the four pixels nearest the 
remnant, there is some random noise in the background 
value. We believe this disadvantage is far outweighed 
by the difficulty in doing any kind of more complex fit 
given the spatial intensity variations of the background. 
In addition, the column based interpolation also tends to 
remove any residual striping in the LL cubes that results 
from the effects of rogue detector pixels. 

A plot of the background spectrum is shown in the 
bottom panel of Figure [3l Subtracting the background 
eliminates the contributions of [S III], [Si II] and PAHs 
entirely from the SNR spectrum showing that emission 
from these species is not related to the remnant. All 
of these emission features are typical for the ISM irra- 
diated by the average interstellar radiation field (ISRF). 
The local background spectrum is slightly negative at 
short wavelengths. This is most likely due to the gradi- 
ent in the zodiacal and Galactic foregrounds between the 
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Fig. 2. — Results of the convolution described in Section l2.3l for a bright point source in the N22 region. The panels show the point source 
at (1) 23.0 fim, (2) 30.0 fim, and (3) 36.7 /im. The bottom panels (lb), (2b) and (3b), show the point source before the convolution. The 
top panels, labeled (la), (2a) and (3a) show the source after convolution. The image stretch is logarithmic to highlight the faint diffraction 
rings around the point source. After the convolution, the diffraction rings are aligned at all wavelengths, showing the effectiveness of our 
sTinyTim derived convolution kernels. 



location of the dedicated "off" position and the remnant. 
Subtracting the local background, however, corrects for 
this over-subtraction in the spectrum of the remnant. In 
the continuum, background emission accounts for 40 to 
60% of the emission. For the emission lines, the back- 
ground subtraction removes a few percent of the inte- 
grated strength of [Ne III] and [O IV]. 

In general, the local emission in each wavelength slice 
has a very similar spatial structure, except at the wave- 
lengths of strong emission lines. In the strong emission 
line from [O IV] and/or [Fe II] at 25.8 fj,m, a diffuse 
halo is evident around the supernova remnant. This 
ionized gas cannot be from material that is collision- 
ally interacting with the remnant since it is well out- 
side the extent of the forward shock. A similar structure 
is seen at the wavelengths of the other spectral lines. 
The ionized material surrounding the r emnant has been 
obser ved previous ly in optical s tudies (iTuohv fc Dopital 
Il983t iBlair et all |2000; Finkcls tein et all l2006f) ! The 
source of the ioniz ation for this surroundin g material has 
been discussed bv lTuohv fc Dopital ()1983f ). It is possible 
that the material was ionized by the supernova shock 
breakout or is currently being ionized by emission from 
the radiative shocks in dense knots of ejecta in E 0102. 
The surroundings of the remnant will be discussed fur- 
ther in Section B31 

3. ANALYSIS 

3.1. The Mid-IR Spectrum 

3.1.1. Sources of Dust Continuum Emission 

We now consider the possible sources for dust 
continuum emission from SNR E 0102: circumstel- 
lar/interstellar dust in the forward shocked region, newly 
formed dust in the reverse shocked ejecta or newly 
formed dust in the unshocked ejecta. The contribution 
of each of these regions to the spectrum depends on the 
mass, temperature and composition of the dust present. 



The dominant heating process for dust in the shocked 
ejecta and shocked CSM/ISM will be collis ions with elec- 
trons (Dwek& Werner 198lt lArendt et a l. 1999). The 
equilibrium temperature of a dust grain that is collision- 
ally heated in a plasma can be found by equating the col- 
lisional heating rate of the grain with its radiative cooling 
rate. The collisional heating rate for a dust grain with 
size a by particles with mass m , density n and tempera- 
ture T is given bv lDwekl (|1987D : 

H{a, T, n) = {32/TTm)^^^na'^n{kT f^^h{a, T). (2) 

Here h{a,T) is a dimensionless parameter which de- 
scribes the grain "heating efSciency" and depends on the 
fractional amount of the incident energy of the p article 
that is deposited in the grain by the collision. iDwekl 
(|1987[ ) develop an approximation for h{a,T) based on 
a comparison of the effective thickness of the grain 4a/3 
and the experimentally measured stopping distance of an 
electron or ion in a solid. We use their formulae in this 
analysis. The radiative cooling rate of the grain is given 
by: 

Lgr{a,Td) = Ana^aT^ < Q{a,Td) >, (3) 

where < Q{a,Td) > is the Planck-averaged value of the 
absorption coefficient for a dust grain of the given size 
and temperature. Equating these heating and cooling 
rates yields the equilibrium temperature of the grain. 

In the remnant, shocked CSM/ISM dust would be lo- 
cated in a shell with an outer radius of 22" and an inner 
radius at the contact discontinuity, the exact position 
of which has not been determined. X-ray spectroscopy 
of the blast wave region shows post-shock material with 
abundances consistent with the SMC, electron tempera- 
tures around 1 keV, and ionization timesc ales r = Jiet 
3 X 10^° s cin -3 (Hayashi et al. 1994; Hug hes et aLll2000l : 
iSasaki et al.l [200Q) . Assuming a timescale of ~ 1000 
years (the age of the remnant) gives electron densities in 
the forward shock region of rig ~ 1 cm"'^ . CoUisionally 
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Fig. 3. — This figure shows (a) the spectrum of SNR E 0102 after the subtraction of the local background and (b) the local background 
spectrum. The spectrum of the supernova remnant shows dust continuum and emission lines from neon and oxygen. The background 
spectrum shows lines from sulfur, silicon, neon and oxygen as well as emission from polycyclic aromatic hydrocarbons around 8 /jm and 
dust continuum. The contributions from silicon, sulfur and PAHs are entirely removed from the remnant's spectrum after background 
subtraction. The background spectrum is slightly negative at short wavelengths due to oversubtraction of the Galactic and zodiacal light 
foregrounds. The background subtraction corrects this oversubtraction in the spectrum of the remnant. 



heated dust grains in this region would have equilibrium 
temperatures of between 30 and 70 K for silicate and car- 
bonaceous grains with sizes between 0.001 and 1.0 /J,m 
(assuming the optical c onstants for in terstellar graphite 
and silicate grains from iLaor fc Drai nc 1993). Figure H 
shows the equilibrium temperature of dust grains in the 
forward shocked plasma as a function of grain size. The 
collision rates in this plasma, as well as in the reverse 
shock, are low enough that the grains are not heated 
above their sublimation temperature and destruction of 
the grains via evaporation. 

Shocked ejecta dust would be located in a shell with 
an outer radius at the contact discontinuity and an in- 
ner radius at 15" . The heating mechanism of the dust 
in the shocked ejecta depends on its density structure. 
When the reverse shock encounters more diffuse parts of 
the ejecta its passage creates hot x-ray emitting plasma. 
X-ray spectroscopy of the shocked ejecta shows a plasma 
composed of oxygen, neon, magnesium and silicon with 
electron temperatures around 0.4 keV (5 x 10® K) and 
ionization timescales of riet ^ 10^^ s cm~'^. Assuming 
t ~ 1000 years, as before, gives densities in the x-ray 
emitting reverse shocked ejecta of ^ 20 cm^'^. Fig- 
ure |4] also shows the equilibrium temperature of collision- 
ally heated amorphous carbon and forsterite (Mg2Si04) 
grains as a function of size given th ese plasma param- 
eters, using optical constants fr om iRouleau fc Martini 
(|1991f ) for amorphous carbon and lJager et al.l ( 20031 ) for 



Amorphous Carbon 

Forsterite 
- - - Graphite 
Sihcate 




0.01 0.10 1.00 

Grain Size (/im) 

Fig. 4. — The equilibrium temperature of coUisionally heated 
grains of different compositions in the forward and reverse shocked 
x-ray emitting plasma and the hot, post-shock region of a dense 
knot. We plot only amorphous carbon and forsterite in the reverse 
sho cked ejecta based on the results of our dust model fits in Sec- 
tion l3.4l We assume that dust grains in the shocked CSM/ISM will 
have the optical properties characteristic of interstellar graphite 
and silicate from .Laor fc Draine. (il993i) . 
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forsterite and assuming the heat ing efficienci es for car- 
bonaceous and silicate dust from lDwe3 (|1987t ). 

The interaction of the reverse shock with material in 
dense clumps of ejecta drives radiative shocks that pro- 
duce optical and infrared line emission, discussed further 
in Section [3. 1.2( and heat the dust grains by a combina- 
tion o f radiative and coUisional processes. lArendt et al.l 
(|1999r ) compared the coUisional and radiative heating for 
dust in the hot, post-shock region of a knot in Cas A 
assuming the density, temperat ure and cooling function 
predicted by the models of Sutherland fc Dopital ()1995D 
for shock velocities between 150 — 200 km s~^ — the same 
models which are found to best reproduce the optical 
and UV emission from E 0102. They found that the 
coUisional heating exceeds radiative heating under these 
conditions by factors of a few hundred. Thus, for dust in 
the immediate post-shock part of dense knots, we can 
ignore radiative heating. Figure [4] shows the equilib- 
rium temperatures for coUisionally heated dust grains in 
this region, using Tg = lO'' '''^ K, Ue = 400 cm"^ and 
the gas cooling function F 10~^^'^ ergs cm'^ s~^ from 
[Sutherland fc Dopital (|1995D . On the other hand, in the 
cooled post-shock gas of the dense knots and the in cold 
pre-shock gas, the dust is likely he ated by the radiation 
field generated by the shock front. iBouchet et al.l ()2006[ ) 
calculate the radiative heating for dust in similar con- 
ditions in SN 1987A and find that it can reach temper- 
atures of ~ 125 K. Given the angular resolution of our 
observations we will not be able to distinguish regions 
where dust in knots is coUisionally or radiatively heated, 
but we note that the equilibrium dust temperatures in 
these circumstances are all between ^ 100 — 200 K. 

If the dust grains in these various coUisionally heated 
regions are small enough, a collision with a single elec- 
tron may be sufficient to raise the grain temperature well 
above the equilibrium value for a short time. After this 
temperature spike the grain can cool to below the equilib- 
rium temperature, resulting in a large range of tempera- 
tures for a given grain size. An important consideration 
when dealing with this stochastically heated dust is that 
only some fraction of the grains will be warm at a given 
time, so our mid-IR observations may miss a substantial 
fraction of the dust. Grains will be stochastically heated 
in this manner when two conditions are met: 1) the coUi- 
sional rate is slow compared with the cooling rate of the 
grain and 2) the amount of energy deposited by a single 
electron is large compared with the internal energy of the 
grain. 

For the plasma conditions in the forward and reverse 
shocks. Figure [5] shows the ranges of grain size for which 
we expect the dust to be in equilibrium with the plasma 
based on the ratios between the energy deposited by one 
electron and the grain's total internal energy (AE/U) 
and the ratio between the cooling time and collision 
time {tcooI / Tcoii) ■ For simplicity we show only one grain 
species for the forward and reverse shocks: graphite in 
the forward shock and amorphous carbon in the reverse 
shock. The size at which equilibrium is reached is a fac- 
tor of ~ 2 smaller for silicate grains in the forward shock 
and a factor of ~ 1.5 smaller for forsterite grains in the 
reverse shock. For this calculatio n we have used the heat 
capacity of graphite grains fr o m |Dwekl ([.1986) and sili- 
cate from iDraine fc AndersonI (|1985D . both of which are 
based on fits to experimental data, and the same optical 
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Fig. 5. — This plot illustrates the range of dust grain sizes that 
will be in equilibrium in the forward and reverse shocks. The dotted 
lines are for amorphous carbon in the reverse shock and the dashed 
lines are for graphite in the forward shock, assuming the densities 
and temperatures discussed in the text. Two ratios are shown 
on the plot; the ratio between the cooling time and the collision 
time (t^ooI / '''coll) S'lid the ratio between the energy deposited by 
a collision with one electron and the internal energy of the grain 
(AE/U). When the ratio of the collision time to the cooling time 
is large (i.e. collisions happen faster than the dust grain can cool) 
and the ratio of the energy deposited by one electron to the total 
internal energy of the grain is small (i.e. the grain temperature 
does not change appreciably from one collision) then the grains 
can be said to be in equilibrium. The grain size for which these 
conditions are met are > 0.01 fim for amorphous carbon in the 
reverse shock and > 0.02 fim for graphite in the forward shock. 
This size is a factor of ~ 1.5 smaller for forsterite in the reverse 
shock and ~ 2 smaller for silicate dust in the forward shock. 

constants used in Figure [H The graphite heat capacity 
is used for amorphous carbon and the silicate heat ca- 
pacity is used for forsterite. For the estimates presented 
here, this level of accuracy should be sufficient, though 
the heat capacity of a dust grain does depend on its com- 
position, shape and "Uuffiness" among other factors. 

The unshocked ejecta are located inside the reverse 
shock radius of 15" . The density, temperature and struc- 
ture of the unshocked ejecta are unobserved. We can 
place limits on the temperature of the gas and dust 
by considering the processes responsible for heating and 
cooling the ejecta. In the initial phases of the remnant's 
expansion the ejecta are heated by the decay of radioac- 
tive isotopes produced in the explosion. As the remnant 
expands, the ejecta cool adiabatically, reaching a temper- 
ature cool enough to allow dust c ondensation after about 
400 to 800 da ys post-explosion (|Todini fc Ferrarall200ll : 
iNozawa et al.l [2003). Assum ing adiabatic coolin g with 
an index 7 = 1.25 as used by iKozasa et al.l (|1989l ) and a 
temperature of ^ 1000 K 400 days after the supernova, 
the current gas temperature would be less than the tem- 
perature of the cosmic microwave background. Before 
reaching that temperature, of course, the dust and gas 
would come into equilibrium with the local radiation field 
which has contributions from both the general SMC ra- 
diation field and the emission from E 0102's shock fronts. 
The equilibrium temperature of big grains in the inter- 
steUar radiation field of the SMC is 20 K (jBot et al.l 
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[200I iLerov etsIMOTh . so we consider this a lower limit 
on the temperature of the unshocked ejecta dust. 

3.1.2. Oxygen and Neon Fine- Structure Emission 

Lines 

We observe emission from fine-structure transitions of 
[Ne II], [Ne III], [Ne V] and [O IV] in the mid-IR spectrum 
of E 0102. At the spectral resolution of the LLl order, 
the [O IV] line overlaps with [Fe II] at 25.988 ^m. How- 
ever, given the predominance of oxygen in the supernova 
remnant as observed in the optical, x-ray and UV and 
the dearth of iron emission observed at any wavelength, 
we consider it very likely that the 25.8 fim emission line 
is due to [O IV] . Table [1] lists the detected emission lines 
and their strengths. Examination of the spatially re- 
solved velocity information shows a range of velocities in 
agreement with previous spectroscopic measurements of 
optical line emission. 

The contribution of [O IV] to the 24 iim flux was one 
of the major uncertainties for IStanimirovic et ahl (j2005l ) 
in their analysis of the mass of newly formed dust in the 
remnant. We find that [O IV] is responsible for ~ 15% of 
the 24 fim flux from the remnant, significantly less than 
the 60% they were forced to assume due to the lack of 
spectroscopic information. 

The line emission from oxygen-rich supernova rem- 
nants like E 0102 is thought to arise from slower, ra- 
diative shocks created when the reverse shock encoun- 
ters dense clumps of ejecta. Due to the high abun- 
dance of oxygen and other heavy elements in the nuclear- 
processed stellar ejecta, the cooling time behind these 
shocks is shorter than the recombination time, produc- 
ing a plasma that is much cooler than its ionization state 
would indicate. The relative ionization levels in the post- 
shock ejecta depend upon the details of the shock front 
and the ionization state of the material prior to encoun - 
tering the reverse shock. [Sutherland fc Dopital (|1995l ) 
found that the optical emission spectrum of E 0102 could 
be fit relatively well with a model for the reverse shock 
entering a dense ejecta knot with velocities between 150 
and 200 km s^^ and creating a precursor photoionization 
front that determines the initial ionization state of the 
pre-shock gas. Using the observed mid-IR emission line 
strengths to determine physical conditions in the ejecta 
would require shock modeling outside the scope of this 
paper. 

3.2. Spatial Decomposition of the Mid-IR Emission 

At the resolution and signal-to-noisc of our observa- 
tions it is difficult to spatially separate the emission from 
circumstellar/intcrstellar dust in the forward shock from 
the emission from newly formed dust heated by the re- 
verse shock. The outer radius o f the forward shock, as 
measured by iGaetz et al.l ()2000t) is 22" while the in- 
ner radius of the reverse shock is ^ 15". In between, 
the two plasmas are separated by a contact discontinu- 
ity, the location of which has not been measured. The 
resolution of the IRS instrument in the LLl order ranges 
between 6 and 9" and the pixel scale is ^ 5" . The angular 
resolution inhibits our ability to separate newly formed 
dust from dust heated by the forward shock, without 
more detailed modeling. Measuring the amount of for- 
ward shocked dust is also an interesting quantity on its 
own (Section [43]). 



The close correspondence between the 0.3 — 10 keV 
Chandra map and the 2 4 /um dust map, as noted by 
IStanimirovic et ahl (|2005D . suggests that the dust emis- 
sion is mainly from the reverse shocked material and is 
thus newly formed in the remnant. In this section, we 
quantify the fraction of the dust emission at each wave- 
length in our cube that can be attributed to dust in the 
ejecta of the supernova by decomposing the remnant spa- 
tially using x-ray, radio and optical observations as tem- 
plates for the sources of emission in the remnant. As 
discussed in the previous section, there are three mecha- 
nisms for which we might expect to detect mid-IR dust 
emission: CSM/ISM dust in the forward shocked region, 
newly formed dust in the x-ray emitting regions of re- 
verse shocked gas and newly formed dust in the dense 
knots of ejecta. Thi s approach is similar to that used by 
lArendt et al.l (|1992D in studying dust emission from the 
Cygnus Loop. We will first discuss the template images 
we use to represent these three emission mechanisms. 

Emission from the forward shocked material can be 
traced by radio continuum emission. The radio contin- 
uum is synchrotron emission from electrons interacting 
with the compressed CSM/ISM magnetic field in the re- 
gion between_tlie blast w ave and the contact discontinu- 
ity. lAmv fc Ba ll (1993) observed the remnant at 6 cm 
using the Australia Telescope Compact Array. They see 
a shell structure which fills the region between the faint 
outer boundary of the blast wave as traced by x-ray emis- 
sion and the outer boundary of the reverse shock. We use 
the radio map as a tracer of the forward shocked circum- 
stellar and interstellar material. 

We have use d the 0.3 — 10 keV Chandra x-ray map 
of IGaetz et"all ()2000l ) to trace the x-ray emitting re- 
verse shocked material. X-ray spectroscopy of E 0102 
has demonstrated that the bright x-ray ring contains 
material composed primarily of oxygen and neon with 
smaller contributions from magnesium and silicon, con- 
siste nt with the nucleosynthetic products of a massive 
star (jFlanagan et al.|[20M : ISasaki et al.ll2006D . Finally, 
to trace the spatial distribution of dust in dense knots of 
reverse shocked material we use the optical [O III] 5007 
A emi ssion from HST ACS imaging bv lFinkelstein et al.l 
(|2006f l that arises from the aforementioned radiative 
shocks. 

The templates were convolved to the resolution of the 
longest wavelength point spread function of the IRS cube 
and then resampled to the same pixel grid. For the 
[O III] template a few additional reduction steps were 
necessary prior to the convolving and alignment. The 
processed line (F475W) and continuum (F550M) images 
from HST were provided to us by S. Finkelstein. The 
wide F475W filter was used to ensure detection of the 
high velocity emission from the remnant. Relative to 
the point sources in the image, the remnant emission is 
quite faint. Thus, to properly construct a template it 
was necessary to remove the point sources from the im- 
age by subtracting the continuum image from the line 
image. This was done by first cross-convolving th e two 
images with theoretical ACS PSFs from TinyTim (jKristl 
hggg) for the given wavebands, scaling, and subtracting. 
For bright point sources, the PSF subtraction often left 
large residuals, which were removed with a combination 
of median filtering and masking by hand for the bright- 
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TABLE 1 

Emission Lines in the Mid-IR Spectrum of E 0102 



Line 


Measured Wavelength 


Central Intensity 


Inte 


grated Strength 




(/xm) 


(MJy sr-i) 


(10- 


7 W m-2 sr-i) 


[Ne II] 


12.8053 ± 0.0005 


6.287 ± 0.035 




2.585 ± 0.015 


[Ne III] 


15.5514 ± 0.0010 


2.295 ± 0.019 




0.917 ± 0.022 


[Ne V] 


24.3460 ± 0.0053 


0.566 ± 0.015 




0.234 ± 0.017 


[O IV] 


25.9022 ± 0.0010 


4.662 ± 0.019 




2.132 ± 0.014 



Note. — These measured values are from Gaussian fits to the IRS spectra. 



est, saturated stars. The three templates are shown in 
Figure[n]after convolution and resampling to the LL pixel 
grid. 

The decomposition at each wavelength was done with 
a least squares fit that constrained the coefficient of each 
template image to be greater than or equal to zero. The 
fit returns the best coefficients of the linear combination 
of the three templates that match the observed surface 
brightness distribution. Examples of the decomposition 
are shown in Figures [7] and [5] for two images which show 
20.5 — 22.3 fj,m dust continuum and the [O IV] emission 
line at 25.9 /xm, respectively. To obtain the images used 
for these plots we have binned the cube in wavelength to 
increase the signal-to- noise of the image. Figure [9] shows 
the results of the decomposition for all wavelengths. The 
error bars on the plot include covariance between the 
three templates. One can interpret each of these re- 
sults as the spectrum of the emission from the region 
described by each template: the x-ray emitting reverse- 
shocked ejecta, the [O III] emitting reverse shocked ejecta 
in dense knots and the forward shocked CS]\/I/IS]VI ma- 
terial. 

It is clear from the residuals shown in Figures [7] and [8] 
that the correspondence between the templates and the 
observations, while good, misses some features of the in- 
frared emission. There is a persistently bright region in 
the lower left of the images which has no counterpart 
in the combined template. At the wavelength of strong 
emission lines the combined template subtraction leaves 
a ring of emission which most likely represents an un- 
dersubtraction of the background at these wavelengths. 
Because of these issues, when determining the best fit 
dust mass we consider the most reliable answer to be 
from the fit to the total remnant spectrum, rather than 
the combination of the spectra associated with x-ray and 
[O III] emitting ejecta. 

Despite these uncertainties there are some strong con- 
clusions that can be drawn from the decomposition re- 
sults. One results that is a good test for the technique 
is that all of the line emission is found to come from 
the regions traced by optical [O III] emission, as ex- 
pected. Other results of the decomposition are: 1) that 
the spectrum associated with the x-ray emitting ejecta 
shows dust continuum but no emission lines, 2) there is 
some dust continuum in the region traced by [O III] emis- 
sion, which seems to peak at a wavelength around 22 fim 
and 3) the forward shock spectrum also has a small dust 
emission component that peaks around 20 /xni. We dis- 
cuss the implications of the forward shocked dust further 
in Section B31 



3.3. Modeling the Dust Emission 
3.3.1. Constructing the Dust Model 

Having established that the mid-IR dust continuum 
emission from E 0102 mainly originates in the reverse 
shocked ejecta of the remnant, we now proceed to model 
the dust emission to attempt to place limits on the total 
amount of dust formed in the CCSN. Because we only 
see emission from material that has been reheated by 
the reverse shock, our model for the dust must include 
only the types of dust that can form out of the nucle- 
osynthetic materials present in the layers of ejecta t hat 
are shocked (iKoza sa et"allll989l: iNozawa et al.ll2003 V A 
crucial observation is that the ejecta in SNR E 0102 
do not appear to be mixed. Observati o ns at x-ray 
(iFlanagan et all 12004 ISasaki e tlHI 120011: [Caetz et all 
[2000; Rasmuss en et al.l 2001: H avashi et al.l 119941) . uF 
traviolet (Sasaki et al. 2006), optical fFinkclstein et al.l 
[2006; Tuohy & Dopita 1983; Blair et al. 2000, 198^) and 
now infrared only detect emission from oxygen, neon, 
carbon, magnesium and silicon. This is in contrast 
to the similar 0-rich remnant Gas A, which shows 
emiss ion from sulfur, argon, ca lcium an d iron a s well 
(Hwa ng et aLll2004l : iFesen et al.l [2006: R ho et all [20081. 
Either the nucleosynthetic layers of the progenitor star 
appear to have undergone some macroscopic mixing in 
Gas A or the reverse shock has encountered the oxygen- 
burning layers and the carbon-burning layers, for in- 
stance, iri_difierent parts of the supernova, simultane- 
ously ([Ennis et al.l I2006D . The same is not true of E 
0102 where we see only nucleosynthetic products from 
the carbon-burning layers. This is most likely a conse- 
quen ce of the ej e cta be ing stratified and unmixed. 

Noz awa et al.l ([20031 ) model the condensation of dust 
in the unmixed ejecta of core-collapse supernovae. The 
types of dust they predict are amorphous carbon in the 
He- rich layers; AI2O3, Mg2Si04 (forsterite) and IVIgO in 
the 0-Mg-Si layer; MgSiOs and Si02 in the 0-Si-]VIg 
layer; and silicon and iron rich species in the deeper nu- 
cleosynthetic layers. To construct our model for the dust 
in E 0102 we first consider how deeply into the ejecta 
the reverse shock has propagated, and thus what species 
of dust we should include in our model. Measurements 
of the abundances of oxygen, neon, magnesium and sil- 
icon in the x-ray spectra of the shocked ejecta demon- 
strate t hat magnesium is 2 times more abundant than 
silicon ([Flanagan et al.l [20041 ). The relative amounts of 
these two elements locate the reverse shock in the 0-Mg- 
Si layer. For this reason, in our dust model we assume the 
primary species are amorphous carbon, AI2O3, forsterite 
and MgO as predicted for the 0-M g-Si layer and the He- 
rich layer bv iNozawa et al.l ([2003f ). 
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Fig. 6. — Convolv ed and al igned x-ray, optical and [O III] templates used in decomposing the mid-IR emission. Data are from lGaetz et aP 
(|2003) in the x-ray. lAmv fc BalL (il993) ') in the radio and iFinkelstein et al.i (|2006|) in the optical. 
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Fig. 7. — Results of the spatial decomposition for a stack of continuum images between 20.5 and 22.3 /im. We show the stacked image 
in order to increase the signal-to-noise ratio for display purposes. Panel a) shows the stack, Panel b) shows the best fit combination of the 
three templates and Panel c) shows the residuals of the fit. Note that the residual image has a different grey scale than the images in the 
Panels a) and b). While the linear combination of the templates reproduces most of the features of the mid-IR images, the residuals do 
have some persistent features. 




Fig. 8. — Same as in Figure [7] but for the [O IV] emission line at 25.9 fim. Note that the grey scale for the residual image has a 
different scale than the other images. The ring of emission visible in the residual image is most likely due to a slight undersubtraction of 
the background in the vicinity of a bright emission line. 
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Fig. 9. — Results of the spatial decomposition of the mid-IR spectrum. Each panel shows the spectrum of the emission associated with 
each tracer. Error bars include the covariance between the tracers. Spikes in the x-ray and radio results in the vicinity of spectral lines 
arc artifacts of the decomposition technique. The dotted line overplotted in each panel shows the total remnant spectrum. The [O III] 
tcmijlatc traces the emission from dense knots of reverse shocked ejecta, the x-ray template traces shocked ejecta and the radio template 
traces shocked CSM/ISM. 
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The importance of this approach to the determina- 
tion of the dust mass in E 0102 arises from trying to 
find a physically motivated and yet unique model for 
the dust composition. Given enough dust species out 
of the complicated mineralogy of dust condensation, one 
could imagine any number of degenerate fits to the mid- 
IR spectrum. By restricting the dust model to amor- 
phous carbon, AI2O3, MgO and forsterite we hope to 
obtain physically meaningful values for the dust mass of 
each species. As we will discuss further below, even with 
these constraints the fit is not entirely unique, demon- 
strating the challenges of dealing with dust continuum 
emission. 

3.3.2. Upper Limit at lOjim 

An important constraint on our determination of the 
dust mass from the mid-IR spectrum is the upper limit 
for the surface brightness of the supernova remnant at 70 
^m from the S'^MC MIPS observations. Contamination 
by the local foreground is higher at longer wavelengths 
because of the increasing contribution from cooler dust 
and the lower resolution at 70 /itm means it is more 
difficult to separate the remnant from the shoulder of 
N 76. For these reasons determining the upper limit is 
not straightforward. There is no obvious emission visible 
from the supernova remnant in the MIPS 70 /im image. 
To place a limit on the remnant's surface brightness we 
subtract a local background using the same technique 
outlined in Section 12.21 with a mask region 40" in ra- 
dius to avoid contamination from the remnant, since the 
angular resolution at 70 /im is 18". The background 
subtraction leaves emission consistent with zero, with an 
error of 1.0 MJy sr~^ in the mean. We conservatively 
adopt a 2-tT upper hmit of 2.0 MJy sr~^ at 70 /zm in the 
dust fitting to be conservative. 

3.4. Dust Model Fit Results 

We performed a non-linear least squares fit to deter- 
mine the best fit dust model including the upper limit at 
70 /J,m. We first fit to the total spectrum of the remnant 
extracted over a 22" radius region from the background 
subtracted cube. Our dust model involves the four grain 
species discussed above with a fixed size of 0.1 /^m. For 
dust grains in the Rayleigh limit the dust mass is inde- 
pendent of the grain size, so fixing the size at 0.1 /xm 
does not affect the fit. The parameters of the model are 
the mass of dust in each species and its temperature. We 
also performed the fit on the total remnant spectrum mi- 
nus the spectrum associated with the decomposed radio 
emission shown in Figure [9] The results from these two 
fits are very similar, as shown in Table [2l Removing the 
spectrum associated with the shocked CSM/ISM slightly 
increases the steepness of the long wavelength continuum 
leading to more cool dust in the fit. 

The spectrum of the remnant is well fit by a combi- 
nation of forsterite (Mg2Si04) and amorphous carbon 
dust as shown in Figure [TOl Neither AI2O3 nor MgO 
play an important role in the best fit model. Upper 
limits on the masses of AI2O3 and MgO are listed in 
Table [2] at a temperature of 70 K. At higher tempera- 
tures the limits on the masses are stricter. The best fit 
forsterite component has a temperature of 145 K while 
the amorphous carbon component, which contains most 
of the dust mass we detect, has a temperature of 70 K. 



However, the fit results for the cool dust component are 
not unique. A fit of similar quality using two tempera- 
ture components of forsterite is listed in Table [5] as well. 
The continuum shape around 20 /im is not well repro- 
duced without forsterite in the dust model, however the 
longer wavelength continuum shape does not distinguish 
strongly between forsterite and amorphous carbon. In 
the following we proceed by assuming that the dust is 
amorphous carbon because 1) it provides a more con- 
servative estimate of the dust mass (i.e. a lower limit) 
and 2) theoretical predictions suggest that there should 
be a substantial amount of amorphous carbon which, if 
present, should all be in the reverse shocked outer layers 
and visible to us in the mid-IR. 

The 70 /im upper limit does not affect the amount 
of forsterite we measure. On the other hand, the limit 
does affect both the mass and temperature of the amor- 
phous carbon dust. Decreasing the limit tends to make 
the amorphous carbon dust warmer and less massive, in- 
creasing it makes the amorphous carbon cooler and more 
massive. We note that for limits below about 1 MJy 
sr~^ the fit no longer converges properly, because the 
limit is inconsistent with the slope of the data past 30 
/tm. We predict that with a small increase in signal-to- 
noise and/or angular resolution the remnant should be 
detectable at 70 /im. 

We also perform individual fits to the spectra associ- 
ated with the x-ray, [O III] and radio templates in or- 
der to learn about the dust in the x-ray emitting reverse 
shocked ejecta, the dense knots of reverse shocked ejecta 
and the CSM/ISM dust in the forward shock. Figure [TT] 
shows the dust continuum associated with the radio and 
[O III] templates and the fit listed in Tabled The con- 
tinuum shape and signal-to-noise in the spectrum asso- 
ciated with the radio template do not strongly constrain 
the species of dust present in the CSM/ISM material. We 
show a representative fit with interstellar graphite from 
iLaor fc Drail^ (|1993f ). The fits imply - 10"^ - 10"^ M© 
of dust with a temperature around 150 — 180 K. We dis- 
cuss this component further in Section|4?5l The dust con- 
tinuum associated with the [O III] emitting knots has a 
better defined shape that is well reproduced by forsterite 
but not amorphous carbon. Fits to that component yield 
a dust mass of '--^ 2 x 10^^ M© of forsterite at 180 K. 

Using a simplified version of the model which involves 
only forsterite and amorphous carbon, we performed the 
fit at each pixel in the map in order to learn about the 
spatial distribution of the dust species. These results are 
shown in Figures fT2l and [T3l The forsterite to amorphous 
carbon ratio tracks the intensity of the x-ray ring of the 
reverse shock. This resemblance may be an indication 
that forsterite is more abundant in the deeper layers of 
the ejecta which are now encountering the reverse shock. 

4. DISCUSSION 

4.1. Newly Formed Dust 

We detect 3 x 10~^ M© of 70 K amorphous carbon 
and 2 x 10~^ M© of 145 K forsterite in the ejecta of E 
0102. The equilibrium temperatures of grains in the x- 
ray emitting plasma are between 50 and 100 K for sizes 
between 1.0 and 0.001 /im. The temperature we measure 
for the amorphous carbon component corresponds to the 
equilibrium temperature for grains with sizes of ^ 0.1 
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TABLE 2 
Fit Results 



Spectrum 


ivioaei 


iViass 
(Mo) 




Temperature 
(K) 


Total Spectrum 


AI2O3 


< 1 X 10" 


6 


> 70 




Am. Carbon 


3 X 10" 


3 


70 




MgO 


< 1 X 10" 


5 


> 70 




Mg2Si04 


2 X 10" 


5 


145 


Total Spectrum 




2 X 10~ 


5 


145 




Mg2Si04 


1 X 10" 


2 


54 


Total Spectrum - CSM/ISM 


Am. Carbon 


6 X 10" 


3 


60 




Mg2Si04 


2 X 10" 


5 


145 


X-ray Emitting 


Am. Carbon 


4 X 10" 


3 


60 


Reverse Shocked Ejecta 


Mg2Si04 


8 X 10" 


6 


150 


[O III] Emitting 


Mg2Si04 


2 X 10" 


6 


180 


Reverse Shocked Ejecta 










CSM/ISM 


ISM Graphite 


8 X 10- 


6 


170 



Not e. — Optical constant s for amorphous carb on from iRouleau Sz Martini 
<]199iri and ISM graphite from [Laor fc Draind lll993f) . The sources for the other 
optical constants can be found in Table [3] 




100 



Wavelength {fim) 



Fig. 10. — Dust model fit to the spectrum showing 3 X 10~^ Mq of 70 K amorphous carbon and 2 X 10~° M-q of 145 K forsterite. The 
145 K forsterite component is necessary to reproduce the continuum shape below ~ 30 iim. The 60 K dust component is less constrained 
and could be replaced with forsterite at a comparable temperature with no impact on the quality of the fit. The arrow shows the upper 
limit at 70 /xm. 
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Wavelength {fJ-m) 



Fig. 11. — This figure shows a zoomed in version of the spectra associated with the [O III] and radio templates from the spatial 
decomposition along with a fit to the dust continuum. The fit to the spect rum associated with th e radio template is representative. Since 
we expect it to be interstellar dust we have used interstellar graphite from lLaor fc Draind ||1993IV Silicate also reproduces the continuum 
shape at the same level of quality. For the [O III] knots, however, the shape of the continuum is not well reproduced by amorphous carbon. 
The fit using forsterite is quite good, given the uncertainties involved. 
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Fig. 12. — Dust masses determined at each pixel. The first panel shows the mass of amorphous carbon, the second shows the mass of 
forsterite and the third shows the ratio of forsterite to amorphous carbon. The ratio is highest at the boundary of the reverse shock which 
agrees with the idea that the reverse shock is just beginning to encounter the O-Mg-Si layer where forsterite has been produced. 
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Fig. 13. — Dust temperatures determined at each pixel. The temperature of the amorphous carbon dust is consistent with grains of ~ 0.1 
fim in equihbrium with the x-ray emitting reverse-shocked plasma. The forsterite component has temperatures which are higher than the 
predicted equilibrium in the x-ray emitting plasma. This may be the result of th e forsterite abundance increasing deeper in the ejecta, as 
suggested by the forsterite to amorphous carbon abundance shown in Figure [121 where the recent passage of the reverse shock has not yet 
destroyed smaller grains which can be stochastically heated. 



fim. Figure [5] verifies that grains of this size should be in 
equilibrium. 

The temperature of the forsterite dust component is 
hotter than the equilibrium temperature that could be 
achieved by grains in the x-ray emitting plasma of the 
reverse shock, but is consistent with temperatures in the 
dense [O III] emitting knots. In fitting the decomposed 
spectrum associated with the [O III] emission we do see 
forsterite with temperatures around 180 K, typical for 
grains with sizes around 0.01 fim under the expected 
conditions. If, however, most of the forsterite is located 
in the x-ray emitting plasma, as the results of our spa- 
tial decomposition suggest, the temperature would imply 
small dust grains that are stochastically heated. In the 
following section we show that the sputtering timescale 
in the reverse shocked ejecta for small grains is less than 
the age of the remnant. The presence of small, stochasti- 
cally heated dust grains in the x-ray emitting plasma can 
be explained if the dust has only recently been through 
the reverse shock, such that it has been in the plasma for 
less than the sputtering time. Given the results shown in 
Figure [T^ i.e. that the abundance of forsterite seems to 
increase near the bright, inner boundary of the reverse 
shock, this may very well be the explanation for why we 
see forsterite grains out of equilibrium with the plasma. 

4.2. Dust Sputtering and Grain Size Constraints 

The dust we see in E 0102 at present may not be the 
amount that was initially produced in the supernova or 
what was initially swept up by the blast wave. Dust 
grains in shock waves can be destroyed by sputtering, 
where impacts by hot gas particles erode the grain. 

For the dust in the forward shocked ISM/CSM we can 
estimate the sputtering time with the relation 



' sp 



da/dt 



10' 



fim Vcm^"^/ 



from lDraine &: Salpetej (|1979f) for gas with temperatures 
greater than ~ 10^ K and solar abundances. In this 
regime the majority of the sputtering is done by hydro- 
gen ions and Ue ~ nn, so the abundance of heavier el- 
ements does not significantly affect the calculations and 
we can use these results for the SMC metallicity. Follow- 
ing the discussion in Section [3. 1.11 we estimate densities 
in the x-ray emitting forward shocked material of ~ 1 
cm^^. Thus, the sputtering timescale for dust grains in 
the forward shock is ^ 10^ years for an 0.1 fim grain and 
~ 10'^ years for an 0.001 fim grain. Given the age of the 
remnant, dust grains with sizes below 0.001 fim should 
have been destroyed behind the shock. 

In the reverse shocked gas the sputtering of dust is 
due to collisions with oxygen ions, since oxygen is the 
dominant component of the ejecta. We can use a simple 
scaling argument based on the sputtering time relation 
above to estimate the lifetime of grains in the reverse 
shocked material. The sputtering rate R in (g s~^) is 
given by 

R - crnvY, (5) 

where n and v are the density and velocity of the sput- 
tering particles, a is the grain cross section and Y is the 
sputtering yield per collision in units of mass per colli- 
sion. Therefore, the change of grain radius with time, 
assuming constant density for the grain and a thermal 



velocity for the colliding particles v 
as: 



(kT/m 



,l/2_ 



da „i/2 
— (X ni ' m 
dt 



goes 



(6) 



(4) 



Scaling the sputtering between the forward shock and re- 
verse shock gas involves a number of assumptions. First, 
the density of oxygen ions in the reverse shock is not 
simply related to the electron density, because the ion- 
ization state of the gas is not in equilibrium (in which 
case the number of free electrons donated per oxygen 
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atom would be a predictable function of the tempera- 
ture). Second, the temperature of the oxygen ions is 
not an easily determined quantity either, given the un- 
certain equilibration timescale between the electrons and 
the ions after the shock. As an estimate we can assume 
no ~ ne/3 as a reasonable value for the number of free 
electrons donated by each oxygen atom, and To ^ 7e, 
which should be the case if the equilibration timescale be- 
tween the ions and electrons is less than ~ 1000 years. In 
the case where the ions are hotter t han the electrons, the 
sputtering timescale will decrease. iNozawa et al.l ()2007t ) 
investigated the sputtering of dust grains in the reverse 
shocks of supernova remnants and via fits to experimen- 
tal data for a variety of grain compositions and collider 
properties found that sputtering by oxygen ions has ^ 50 
times higher yield than sputtering by hydrogen ions. Us- 
ing these values (no ~ "e/S, Te ~ To, and Yq ^ 501h) 
we estimate that the sputtering timescale in the reverse 
shocked ejecta, which has an electron density of ~ 20 
cm~^, for a grain of size 0.1 ^vci is ^ 1000 years. Thus, 
in the x-ray emitting, oxygen-rich ejecta, the destruction 
time for grains smaller than 0.1 /im is shorter than the 
age of the remnant and we expect that the grain size 
distribution in the plasma has been altered significantly. 

The amorphous carbon dust, which makes up the ma- 
jority of the mass we detect from the reverse shocked 
ejecta, has a temperature which suggests emission from 
~ 0.1 pjn grains. Based on our estimate of the sputter- 
ing timescale in the reverse shocked plasma, these grains 
would be able to survive the lifetime of the remnant. The 
warmer forsterite dust component is found to be associ- 
ated with the dense [O III] knots, which have a higher 
equilibrium temperature due to their high densities, and 
with the most recently reverse-shocked ejecta, in which 
smaller dust grains may still exist and be stochastically 
heated to higher temperatures. 

4.3. The Possibility of Cold Dust 

Because our observations are only sensitive to warm, 
recently shocked dust, we cannot quantify the total 
amount of newly formed dust produced in the remnant. 
To understand the limits that the 70 /im observation 
places on the mass of cold dust, we construct a model 
for the unshocked dust in the s upernova remnant b ased 
on the condensation models of iNozawa et al.l (|200 3') for 
a 25 Mq progenitor. We list the dust species included 
in the model, the sourc e for their opti c al con stants, and 
the masses predicted bv INozawa et al.l (|2003[ ) in Table H 
The total mass of dust predicted by their models is 0.6 
Mq, minus the amorphous carbon component which we 
argue must primarily be in the outer, shocked layers of 
the ejecta. In Figure [T4l we show the remnant's spectrum, 
the best fit model with amorphous carbon and forsterite 
from Section 13.41 and the spectrum of the cold dust at 
temperatures of 20 (blue) and 35 K (black). It is clear 
that our limit at 70 /im does not preclude the possibility 
of a substantial mass of cold dust in the remnant, since at 
a temperature of 20 K, the entire predicted mass of dust 
can be hidden at 70 /im. Determining the total amount 
of newly formed dust in E 0102 will require longer wave- 
length observations. 

4.4. Implications for Dust Production in CCSN and 
Comparison with Previous Results 



IStanimirovic et all ()2005[ ) found a mass of 8 x 10"'' M© 
of dust in E 0102 based on the flux at 24 /im. Using upper 
limits at 8 and 70 /im, they argued that the dust had to 
have a temperature around 120 K. They assumed that 
60% of the 24 /im emission was from the contribution of 
[O IV], a value we measure to be 15%. These differences 
alone do not ex plain the ^ 4 times more dust that we 
find compared to IStanimirovic et al.l (|2005l ). Instead, the 
greater contribution of cooler dust with a temperature 
of around 70 K accounts for most of the difference. Our 
measurement of the temperature and mass of forsterite 
dust are comparable to their predictions. The spectrum 
past 30 /im, however, is not compatible with only a ~ 120 
K component. 

Similar mid-infrared spectral mapping observations 
have been carried out on Gas A, ano t her yo ung, oxygen- 
rich CCSN remnant, by iRho et al.l (|2008f ). There are 
some important differences between the observations. 
Cas A's relative proximity and brightness allows more 
detailed fits to the spectrum of the remnant. Cas A 
also has emission from S, Si, and Ar in its spectrum in 
addition to O, Ne, and Mg. This coul d be the result 
of ma croscopic mixing, as suggested by iDouvion et al.l 
(|1999D . or of velocity inhomogeneity (Smith et al., in 
prep.), causing different layers of the ejecta to encounter 
the reverse shock simultaneously. In both of these sit- 
uations, the overall composition of the nucleosynthctic 
layers is preserved and the same species of dust should 
have formed in the ejecta as in E 0102. It may not be 
possible, however, to separate different nucleosynthctic 
l ayers spatia l ly, ma king the modeling more complex. 

iRho et al.l (|2008f) find on the order of 0.02 - 0.05 Mq 
of dust in the Cas A remnant, with contributions from a 
variety of dust species. They find three distinct spectra 
in Cas A, a 21 /xm peak spectrum, a weak 21 /im spec- 
trum and a featureless spectrum. The weak 21 /im spec- 
trum has associated strong neon lines, similar to what we 
see i n E0102. The do minant dust components by mass 
that iRho et al] ()2008l ) fit to the weak 21 /im spectrum 
are two temperature components of amorphous carbon 
at - 80 and ~ 200 K totahng ^ 1 - 2 x 10"^ Mq and 
- 0.6 - 1.4 X 10-2 Mq of FeO at - 60 K. The mass of 
amorphous carbon we find is comparable to their results, 
however we have no evidence of FeO in our spectrum and 
it seems unlikely that a substantial amount of FeO could 
be produced, since it requires microscopic mixing of the 
ejecta. Co mparing with the total amount of dust that 
iRho et ahl ()2008 ) find in Cas A, we see similar amounts 
of both amorphous carbon and Mg2Si04, on the order 
of IQ-^ and 10"^ Mq respectively. However, they argue 
for significant amounts of Fe, FeO, FeS and Si02, pro- 
ducing a total dust mass of between 0.02 — 0.05 Mq of 
newly-formed dust — an order-of-magnitude larger than 
what we observe in E0102. The progress of the reverse 
shock into the deeper nucleosynthctic layers of Cas A 
may partially account for the difference in the observed 
dust mass between the two remnants. It is possible that 
similar dust masses could be present but the iron- and 
silicon-rich dust species are in unshocked layers of the 
ejecta of E 0102 and are too cold for us to detect in these 
observations. 

Given the progress of the reverse shock into the ejecta, 
we can claim with good confidence that we should see es- 
sentially all of the amorphous carbon dust that has been 
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TABLE 3 

Cold Dust Model Components 



Species Mass Optical Constants Reference 
(Mq) 

Mg2Si04 0.2 Jager et al. (2003) 

Amorphous Silicon 0.2 Filler (1985) 

Si02 0.07 Philipp (1985) 

FeS 0.07 Begemann et al. (1994) 

MgO 0.06 Roes sler & Huffman (1991) 

Fe 0.02 Lynch & Hunter ( 1991^ 

MgSiOs 0.003 Dorschner et al. (1995) 

AI2O3 0.002 Begemann et al. (1997) 



10.0 : I I I I I I I r 




Wavelength (/xm) 



Fig. 14. — This figure shows predictions for the far-IR spectrum of dust in the r emnant using t he co ld dust model components listed 
in Table [3] The cold dust spectrum contains the total dust mass predicted by the INozawa et al.l II2003I ) models for a 25 Mq progenitor 
excluding amorphous carbon. The dashed line shows this model assuming the cold dust has a temperature of 20 K and the dotted line 
shows the same model with a temperature of 35 K. The black dash-dot and long-dashed lines show the total spectrum of cold dust plus 
the warm dust we detect in the mid-IR. This plot demonstrates that all of the predicted ~ 0.6 Mq of dust can be present in the remnant 
if it has a temperature of 20 K, showing that our current observations alone can not strongly constrain the dust condensation efficiency. 



produced in the mid-IR. This being the case, we find sub- 
stantially less amorphous carbon dust tha n predicted by 
dust c ondensation models. The models of iNozawa et alj 
()2003D predict on the order of 10"^ - 10"^ Mq of amor- 
phous carbon, whereas we only see ~ 3 x 10"'^ M0. 
We estimate that the sputtering time for dust in the 
reverse-shocked ejecta is less than the age of the rem- 
nant, so we may be observing less dust then was initially 
formed. INozawa et al.l (|2007l ) calculate the dust destruc- 
tion via sputtering in the reverse shock for the same mod- 
els quoted above using a simple prescription for the in- 
teraction of the ejecta with the surroimding ISM. They 
find that for an ambient density of 0.1 cm"'^, only 8% 
of the amorphous carbon dust is destroyed for their 25 
M0 progenitor. Most of the mass of carbon dust in that 
model is in grains with sizes greater than 0.05 fim, which 
have longer sputtering lifetimes. We have no constraints 
on the initial grain size, so it is difficult to estimate how 



much of the dust in the remnant has been destroyed up 
to this point. 

Previous observations of newly formed dust in CCSN 
have yielded dust masses in the range of 10"'^ to 10~^ 
Mq, with some controversial measurements claiming 
higher masses. The majority of the measurements come 
from CCSN of Type II or their remnants, which have hy- 
drogen envelopes and red supergiant progenitors. Dust 
formation has been observed in two Type lb events: 
SN 19901 (lElmhamdi et al.ll2C)04D and the peculia r event 
SN 2006ic (jSmith et al.ll2008HNozawa et al.ll2008f ) which 
shows evidence for dust formation at very early times, 
around 50 days after the explosion. The masses deter- 
mined from observations of 2006jc range from a few 10~^ 
to a few 10~^ depending on the ob servation date , wave- 
length and dust model. Recently, INozawa et al.l ()2008D 
have modeled dust condensation in a Type lb event in 
order to assess the effects of a small or non-existent hy- 
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drogen envelope on the dust condensation efRciency. The 
high expansion velocities mean that the density drops 
quickly, leading to dust condensation at earlier times in 
these events. The condensation efficiency is high, how- 
ever, yielding 0.7 Mq of amorphous carbon in their model 
of SN 2006jc, which is in contrast to the relatively small 
amount of dust detected in the near-IR and mid-IR ob- 
servations. Our results for E 0102 lie near the upper 
bounds of the observed dust mass in SN 2006jc. It may 
not be the case, ho wever, that dust is forming in the 
ejecta of SN 2006jc. ISmith et all (|2008f ) argue that the 
dust condensation occurs in the coohng region of the 
shocks created by the interaction of the blast wave with 
an expanding shell of circumstellar material. If this is 
the case, we may not expect to see comparable species 
or amounts of dust in SN 2006jc and E 0102. 

4.5. The Circumstellar Environment of SNR E 0102 

The progenitor of SNR E 0102 was most likely a mas- 
sive star that lost most of its hydrogen/helium envelope 
before exploding as a Type Ib/Ic or Ilb/L supernova 
(jBlair et all 120001 : IChevalieH l2005l ). As such, it would 
have significantly impacted the circumstellar medium in 
the form of winds and mass loss. Here we briefiy dis- 
cuss what we can learn about the CSM/ISM from our 
observations. 

The spectrum associated with the radio continuum 
template, as discussed in Section lX^ shows the dust con- 
tinuum emission from forward shocked CSM/ISM mate- 
rial. Fits to this spectrum with interstellar graphite and 
silicate grains give masses in the range 10~^ — 10~^ M© 
and temperatures of ^ 180 K. If the remnant were ex- 
panding into the average SMC interstellar medium with 
a density of n 0.1 cm~'^, we would expect approxi- 
mately 2 Mq of swept up gas and ~ 3 x lO"*^ Mq of 
swept up dust at this point in ti me (using a dust t o gas 
ratio of 1/700 as determined bv lLerov et al.ll2007| ). Al- 
though the fit dust mass is far less than this expected 
value, the temperature indicates that the dust we see is 
not in coUisional equilibrium and therefore we are only 
detecting the fraction of the dust that is warm at a given 
time (see Figure 2] for the predicted equilibrium temper- 
ature in the forward shock). 

For stochastic heating, the fraction of grains which are 
warm at any time is given approximately by the ratio 
of the collision time to the cooling time / ~ Tcooi/tcoII- 
This ratio is shown as a function of grain size in Fig- 
ure[5l Grains with a > 0.02 fim will not be stochastically 
heated and will have equilibrium temperatures between 
40 and 70 K. In Section 14. 2[ we estimated that grains 
with a < 0.001 /im should be destroyed by sputtering on 
timescales shorter than the age of the remnant. Within 
these limits 0.001 < a < 0.02, we can estimate that 
/ ^ 0.01 — 1, i.e. that we detect between 1 and 100% of 
the dust mass in the forward shocked CSM/ISM. Thus, 
the amount of dust in the forward shocked CSM/ISM 
is essentially unconstrained by our observations, ranging 
from 10-6 - 10-3 Mq. 

Depending on the progenitor star and its winds, we 
may not expect to see any dust in the forward-shocked 
plasma. This would be the case if the blast wave has 
not yet reached the pristine ISM and the circumstellar 
medium has been cleared of du st by the winds of the pro- 
genitor star. iChevalied (l200l argue that E 0102 is still 



interacting with the wind of its progenitor star based 
on the positions of the forward and reverse shocks and 
the lack of limb brightening in the x-ray emission from 
the forward shock. The dust content of the circumstel- 
lar medium will depend on the duration of the wind and 
the previous evolution of the progenitor (for example, 
whether there were outbursts si milar to those observe d 
from the progenitor of SN 2006jc (|Pastorello et al.ll2007[ )) 
among other factors. Further observational constraints 
on the progenitor mass and evolution will be necessary 
to understand whether the dust we see associated with 
the forward shock fits with the overall picture of the evo- 
lution of E 0102. 

5. SUMMARY & CONCLUSIONS 

To summarize, we present mid-IR spectral mapping 
observations of the supernova remnant IE 0102.2— 7219 
in the vicinity of N 76 in the Small Magellanic Cloud. 
The ^ 1000 year old, oxygen-rich remnant is thought 
to be the product of a Type Ib/Ic or Ilb/L supernova. 
Currently the reverse shock, produced by the dynamical 
interaction of the ejecta with the surrounding ISM/CSM, 
is propagating back through the ejecta. Based on spec- 
troscopic observations at x-ray, ultraviolet, optical and 
now infrared wavelengths, the ejecta down to the depth 
of the reverse shock appear to be composed of oxygen, 
neon, magnesium, silicon and carbon. This composition 
indicates that the reverse shock has only penetrated into 
the deeper layers of nucleosynthetic products and that 
no large-scale mixing of deeper layers has occurred. 

The mid-IR spectrum of the remnant shows fine- 
structure emission lines of oxygen and neon on top of 
emission from warm dust. We use observations of the 
remnant at x-ray, radio and optical wavelengths as tem- 
plates to decompose the infrared emission into compo- 
nents associated with the newly-formed dust in the re- 
verse shocked ejecta and ISM/CSM dust in the forward 
shocked region. Our decomposition of the spectrum 
shows (1) dust continuum associated with the x-ray emit- 
ting reverse shocked ejecta, (2) oxygen and neon fine- 
structure emission lines and warm dust continuum asso- 
ciated with the optical [O III] emitting gas in dense knots 
of ejecta where the reverse shock is slow and radiative 
and (3) a small amount of dust emission associated with 
the forward shock. The majority of the dust is associated 
with the x-ray emitting reverse shocked material. 

We construct a model for the newly-formed dust in 
the unmixed ejecta. Using x-ray measurements of the 
abundance of magnesium and silicon we estimate the lo- 
cation of the reverse shock to be in the 0-Mg-Si layer. 
We model the dust emission with amorphous carbon, 
forsterite (Mg2Si04), AI2O3, and MgO as appropriate 
for the composition of the ejecta down to the depth of 
the 0-Mg-Si layer, motivated by the results of dust con- 
densation models. The best fit model contains 3 x IQ-^ 
Mq of amorphous carbon at a temperature of 70 K and 
2 X 10^^ Mq of forsterite with temperature of 145 K. The 
masses of AI2O3 and MgO are limited to be less than 
1 X 10^6 and 1 x 10-^ Mq at a temperature of 70 K, re- 
spectively, since they do not contribute substantially to 
the continuum shape in the mid-IR. The temperature of 
amorphous carbon agrees with expectations for 0.1 /im 
grains in the x-ray emitting plasma behind the reverse 
shock, which agrees with our estimate that smaller grains 



IRS Spectral Mapping of SNR E 0102 



19 



have been destroyed by sputtering. The temperature of 
the forsterite component suggests that the emission is 
from smaller, stochastically heated grains which have re- 
cently been heated by the passage of the reverse shock 
or are residing in dense knots of ejecta material. 

If the reverse shock is presently encountering material 
from the 0-Mg-Si layer of the progenitor we should be 
detecting the majority of the newly- formed amorphous 
carbon dust, since it forms primarily in the outermost 
layers of the ejecta. Comparison with the iNozawa et al.l 
(|2003i) results show substantially less amorphous carbon 
dust in E 0102 than would be predicted based on the pro- 
genitor mass. Our results for the masses of amorphous 
ca rbon and forsterit e are comparable to work on Cas A 
bv lRho et al.' (2008'). 

Our mid-IR observations are not sensitive to cold dust 
present in the remnant. We show that if the dust has 
a temperature of 20 K, ~ 0.6 Mq of dust could be 
present in the remnant and undetectable in our observa- 
tions. Observations at longer wavelengths by Herschel, 
LABOCA on APEX or ALMA will be necessary to de- 
termine the amount of cold dust present in the remnant. 
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